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To  investigate  the  effect  of  energy  on  rock  fragmentation,  initial 
testing  was  carried  out  using  a  drop  test  device.  A  total  of  180  speci¬ 
mens  of  Wausau  quartzite  and  anorthosite  of  3.0  to  3.5  inch  size  were 
fragmented  in  this  device.  -Analysis  of  energy  limitations  of  this  de¬ 
vice  and  other  test  results  led  to  tie  conclusion  that  a  device  supply¬ 
ing  higher  energy  should  also  be  used.  Consequently,  an  impact  pendulum 
test  device  was  designed  and  constructed. 
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20.  ABSTRACT  (Continued) 


An  analysis  of  the  fragment  distribution  resuits  of  the  drop  tests 
showed  that  the  energy  applied  was  approximately  inversely  proportional 
to  the  mean  fragment  si?e.  However,  these  results  were  not  totally  con¬ 
clusive  because  of  the  limited  energy  available  and  the  narrow  band  of 
drop  heights  available. 


A  total  of  270  irregular  3.0  to  3.5  inch  specimens  and  150  disc¬ 
shaped  specimens  of  Wausau  quartzite,  anorthosite,  and  Felch  marble 
were  then  fragmented  with  the  impact  pendulum  device.  -  This  device 
was  instrumented  so  as  to  provide  data  for  calculating  the  specific 
crushing  energy  that  went  into  each  fragmentation  event. ^ Computer 
programs  were  written  to  analyze  sieve  data  and  to  provide  energy 
calculations.  Results  indicated  that  a  power  law  could  be  obtained 
that  provided  the  inverse  relation  between  specific  crushing  energy 
and  mean  fragment  size. 
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FRAGMENTATION  ANALYSIS  -  FUNDAMENTAL  PROCESSES 


INTRODUCTION 

In  the  study  of  rock  fragmentation  there  are  two  fundamental  problem 
areas  which  can  be  investigated.  These  ere  (1)  investigation  of  the 
conditions  under  which  fragmentation  occurs  and  (2)  investigation  of  the 
results  of  fragmenting  processes. 

The  first  of  these  problem  areas  is  where  the  greatest  effort  has 
been  and  is  being  expended.  Most  of  this  work,  however,  should  be  cate¬ 
gorized  as  studies  of  rock  failure  and/or  fracture  rather  than  rock  frag¬ 
mentation.  The  significance  of  rock  fragmentation  is  emphasized  by 
various  efforts  to  relate  failure  mechanisms  to  the  strength  of  rock. 

The  second  problem  area  is  justifiable  from  an  operational  viewpoint. 
Rock  excavation  involves  fragmentation  processes,  such  as  drilling, 
blasting,  and  mechanical  excavation.  Thus,  because  the  nature  of  the 
product  of  a  fragmentation  process  is  a  function  of  the  process  itself, 
and  in  turn  influences  the  selection  of  subsequent  processes  such  an 
further  fragmenting,  handling,  and  transporting  the  rock,  it  is  important 
that  the  nature  of  fragmentation  products  be  known. 

Intuitively,  a  relationship  should  exist  between  a  product  size 
distribution  and  the  energy  input  from  the  nature  of  the  fragmenting 
process  and  the  rock  properties.  Unfortunately,  no  such  relationship 
has  been  defined.  This  then  was  the  purpose  of  this  investigation. 

Specifically,  the  objective  of  this  investigation  was  to  determine 
the  relationship  between  energy  input  and  the  size  distribution  resulting 
from  a  single  elementary  fragmenting  event.  The  actual  investigation  was 
conducted  primarily  to  obtain  the  value  of  the  parameters  that  define 
the  strength  of  the  material  and  the  process  and  to  relate  the  strength 
parameters  with  the  mechanical  properties  of  the  material. 
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TECHNICAL  REPORT  SUMMARY 


The  main  results  of  this  project  were  as  follows:  (1)  The  size  dis¬ 
tributions  of  rock  fragments  resulting  from  a  single  elementary  fragmenta¬ 
tion  event  (comminution  of  a  given  size  rock)  was  obtained.  (2)  This  dis¬ 
tribution  data  was  used  to  determine  the  elements  of  the  breakage  matrix, 
which  relates  the  input  size  distribution  to  the  output  size  distribution. 
(3)  A  power  law  relationship  between  specific  crushing  energy  (or  input 
energy)  and  product  size  was  developed. 

Three  types  of  rocks  were  tested  in  two  types  of  comminution  events — 
drop  tests  and  impact  pendulum  test's.  The  series  c£  drop  tests  were  con¬ 
sidered  to  be  preliminary  tests  to  the  main  series  of  impact  pendulum  tests. 
The  purpose  of  these  dro^  tests  were  twofold— —to  provide  an  alternate 
source  of  comminution  and  to  evaluate  any  possible  effect  of  loading  rate 
on  breakage.  The  loading  rate  had  no  apparent  effect.  Composite  normal 
distributions  were. used  to  fit  the  fragment  distribution  data  from  the 
drop  tests. 

The  impact  pendulum  tests  were  of  the  most  interest.  It  was  found 
possible  to  measure  to  a  reasonable  accuracy  the  velocities  of  the  moving 
parts  before  and  after  impact  so  that  specific  crushing  energy  could  be 
calculated.  A  number  of  different  distributions  were  fitted  to  the  sieve 
data  (resultant  fragment  distributions)  and  it  was  found  that  the  normal 
distribution  provided  the  best  fit. 

From  the  distribution  data,  the  breakage,  matrix,  which  has  elements 
expressing  the  probabilities  of  decreasing  from  one  size  fragment  to  a 
smaller  one,  was  determined.  The  relationship  of  input  (feed)  distribu¬ 
tion  to  output  ^product)  distribution  was  found,  under  certain  assumptions, 
to  be  given  by  p  =  B"  f ,  where  p  is  the  product  vector,  B*  is  the  simpli¬ 
fied  breakage  matrix  and  ?  is  the  feed  vector. 

Relationships  between  specific  crushing  energy  and  average  fragment 
size  were  developed  for  the  impact  pendulum  tests.  For  the  three  rock 
types  tested,  one  such  relationship  was  a  power  law  function  of  the  form 

E/Vq  =  a/yb, 

where 

E/V  **  specific  crushing  energy, 

a ,6  =  constants  related  to  the  breakage  process  and  the  rock  type, 
and  ■  y  =  average  fragment  size. 

An  alternate  relationship  between  specific  crushing  energy  and  average 
fragment  size  is  the  Charles1  law  formulation. 


where 


E/Vq  **  K  [  1  /  X  2 


-  i/xi  n-1L 


E/V  =  specific  crushing  energy, 

K,n  =  constants  related  to  the  breakage  process  and  the  rock  type, 
Xi  *»  initial  specimen  sieve  size, 
and  X2  =  average  product  sieve  size. 

Values  of  a,b  and  K,n  were  determined  from  experimental  data  on  irregular 
and  disc  specimens  of  Wausau  quartzite,  anorthosite,  and  Felch  marble. 
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EXPERIMENTAL  PROGRAM 


Plan  of  Work 


To  simulate  a  single  event  fragmentation  process  that  is  similar 
to  a  real  fragmentation  process  such  as  occurs  with  a  crusher,  it  was 
decided  to  use  random-shaped  specimens  having  a  single  given  initial 
sieve  size.  Disc-shaped  specimens  that  were  relatively  defect-free  were 
also  used  for  greater  homogeneity  of  breakage  with  fewer  specimens. 

Also  the  breakage  results  of  the  disc-shaped  specimens  could  be  compared 
with  the  breakage  results  of  the  random-shaped  specimens.  This  compari¬ 
son  was  essential  because  it  shows  the  two  extremes  of  brittle  fracture 
(see  also  (11)) . 


A  preliminary  experiment  (see  DATA  ANALYSIS  section)  was  designed 
to  determine  the  total  number  of  tests  required  to  achieve  average 
values  representative  of  sample  populations.  Another  test  series  (see 
DATA  ANALYSIS  section)  was  designed  to  study  the  effect  of  specimen  size 
and  shape  for  the  size  range  tested  in  this  test  program. 


Drop  test  equipment  was  initially  used  to  provide  some  of  the  data 
necessary  to  achieve  the  objectives  of  this  test  program.  It  was  known 
that  the  drop  test  data  would  not  account  for  the  energy  dissipated  into 
the  impact  plate,  the  kinetic  energy  utilized  in  scattering  broken  frag¬ 
ments,  heat,  or  acoustic  energy.  Therefore  a  low  velocity  impact  pendu¬ 
lum  was  later  designed  to  improve  on  the  drop  tests. 


For  the  impact  pendulum  tests,  (see  DATA  ANALYSIS  section)  most  of 
the  excessive  energy  available  is  utilized  by  the  moving  parts  of  the 
test  apparatus  and  is  accounted  for  by  (1)  determining  the  velocity  of  the 
impact  pendulum,  first  piston,  specimen,  second  piston,  and  second  pen¬ 
dulum,  and  (2)  accounting  for  heat,  vibration,  acoustic,  and  friction 
energy  by  a  calibration  (no  specimen;  test  procedure. 


Drop  Test  Apparatus 


A  drop  test  apparatus  was  fabricated  and  installed  as  shown  in 
figure  1.  To  minimize  secondary  breakage  (multiple  comminution  of  some 
fragments) ,  the  inside  of  this  chamber  was  lined  with  foam  and  polyeth¬ 
ylene  sheet.  The  specimen  was  placed  on  top  of  trap  doors  in  a  platform 
and  hoisted  to  the  desired  height.  As  soon  as  the  platform  reached  a 
specified  height,  the  specimen  was  discharged  by  a  mechanical  release 
pin  attached  to  the  trapdoor.  Three  drop  heights  of  25,  30,  and  35  ft 
were  selected  to  provide  impact  velocities  of  40.1,  43.9,  and  47.4  ft/sec 
and  impact  energies  of  approximately  62.5,’  75.0,  and  87.5  ft-lb  for  a 
typical  2.5  lb  specimen.  The  broken  material  was  retained  in  a  wooden 
impact  chamber  and  later  was  swept,  bagged,  labelled,  and  sieved  manually. 
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Figure  1.  Drop  test  apparatus 


One  hundred  and  eighty  specimens  (ninety  of  Wausau  quartzite  and 
ninety  of  anorthosite)  were  then  fragmented  in  the  major  test  series 
for  these  drop  tests.  Thirty  specimens  of  each  rock  type  were  dropped 
from  heights  of  25,  30,  and  35  ft.  The  fragments  from  each  specimen 
were  then  sieved  to  determine  composite  size  distribution  of  fragments 
in  each  sieve  class. 


Impact  Pendulum  Test  Apparatus 


A  low  velocity  impact  pendulum  (Fig.  2)  was  constructed  to  provide 
higher  energy  for  rock  fragmentation  and  to  investigate  the  effect  of 
varying  input  energy  on  size  distribution  of  broken  material .  It  mainly 
consisted  of  an  impact  pendulum  (162.5  lb),  first  piston  (38.9  lb), 
second  piston  (41.5  lb),  and  second  pendulum  (269.0  lb).  The  pistons 
were  supported  on  two  sets  of  roller  bearings.  The  parts  of  the  pistons 
were  enclosed  in  the  impact  chamber.  The  test  specimens  were  suspended 
by  a  string  from  the  ceiling  of  this  chamber  to  rest  between  two  piston 
heads  (platens) . 


The  test  procedure  consisted  of  (1)  placing  a  specimen  enclosed  in 
a  plastic  bag  between  the  platens  (the  bag  overlaps  the  platens)  ,  (2) 
releasing  the  impact  pendulum,  (3)  recording  the  horizontal  distance  the 
rebound  pendulum  has  travelled,  (4)  arresting  the  motion  of  both  pendulums 
after  impact  and  travel  measurement  have  taken  place,  (5)  recording  the 
timer  reading. 


Test  Specimens 


Rocks  of  potential  use  in  the  present  investigation  may  be  such  that 
their  constituents  and  properties  vary  considerably  from  one  specimen  to 
another.  To  minimize  such  variations,  three  monomineralic  rocks  -  Wausau 
quartzite,  anorthosite,  and  Felch  marble  were  selected.  The  minerals  that 
these  rocks  are  composed  of  represent  a  large  proportion  of  rock  forming 
minerals  found  in  tunneling  and  other  excavation  projects. 


Irregular  specimens  ranging  in  size  from  3  in  to  3.5  in  were  acquired 
from  muck  piles  of  quarry  blasts.  The  specimens  which  did  not  meet  visual 
inspection  standards  for  uniformity  of  mineral  composition  and  absence  of 
macroscopic  flaws  were  rejected.  Selected  specimens  were  washed,  dried, 
and  weighed. 


In  order  to  determine  the  physical  properties  of  these  rocks,  a 
block  of  each  rock  type  measuring  18  by  18  by  18  in  was  acquired,  fr°in 
which  10  specimens  were  cut.  These  physical  properties  are  presented  in 
Table  1.  Disc  specimens  (2  in  diameter  x  1  in  thick)  were  also  prepared 
from  these  blocks  by  cutting  and  coring.  Again  badly  flawed  specimens 
were  rejected. 
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TABLE  1.  -  Physical  properties  of  Wausau  quartzite, 
anorthosite,  and  Felch  marble 


Rock  type 

Property 

Wausau 

quartzite 

Anorthosite 

Felch 

marble 

Compressive  strength . MN/m2 

. psi 

285.7 

41,437. 

221. 

32,053. 

172.50 

25,019, 

Tensile  strength . MN/m2 

. psi 

5.12 

742. 

8.97 

1301. 

6.78 

983. 

Shore  hardness . 

111.6 

91.7 

65.1 

Pulse  velocity . m/sec 

. ft/sec 

5130. 

16,831. 

6866. 

22,526. 

6686. 

21,936. 

. ft/sec 

4679. 

15,351, 

5,742. 

18,839. 

5,952. 

19,528. 

Torsional  velocity . m/sec 

. ft/sec 

3383. 

11,099. 

3,396. 

11,142. 

3,677. 

12,064. 

Static  Young's  modulus . MN/m2 

. psi 

7.249xl04 

10.51xl06 

4.085xl04 

5.925xl06 

4.318xl04 

6.263xl06 

Dynamic  Young's  modulus*. .MN/m2 

. psi 

5.74xl04 

8.33xl06 

8.971xl04 

13.01xl06 

10.35xl04 

15.01xl06 

Dynamic  shear  modulus** . MN/m2 

. psi 

3.000xl04 

4 . 351xl06 

3.115xl04 

4.518xl06 

3.95xl04 

5.73xl06 

.251 

.310 

.268 

. lb/ft3 

2.649 

165.4 

2.699 

168.5 

2.899 

181.0 

* 

** 

*** 


Calculated  from  bar  velocity 
Calculated  from  torsional  velocity 

Derived  from  ratio  of  longitudinal  pulse  velocity  to  longitudinal 
bar  velocity 
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Figure  2.  Impact  pendulum  apparatus 


Wausau  quartzite  is  a  Precambrian  metamorphic  rock  exposed  in  the 
Rib  Hill  Mountains  southwest  of  Wausau,  Wisconsin,  and  it  is  commercially 
known  as  Wausau  quartzite.  It  consists  of  98.2  percent  quartzite,  0.5 
percent  coarse  amphibole,  0.4  percent  fine  amphibole,  and  0.4  percent 
iron  oxide,  and  0.1  percent  feldspar.  The  impurities  are  scattered 
throughout  the  rock  at  quartz  grain  boundaries.  The  quartz  crystals 
vary  in  diameter  from  3mm  to  8mm. 

Anorthosite  is  a  monomineralic  rock  intrusion  in  the  Duluth  Gabbro 
complex  nearly  5  miles  north  of  Beaver  Bay,  Minn.  It  is  coarse  grained, 
light  green  in  color  and  consists  of  more  than  95  percent  plagioclase. 

Felch  marble  is  a  dolomite  named  for  a  typical  development  at  Felch 
Mountain  in  Dickinson  County,  Michigan.  It  is  mined  from  a  vein  deposit 
intrusion  in  Precambrian  sandstone  by  underground  mining.  The  marble  is 
iti  color*  ecar**.  In  *  and  has  frtmal''  Impurities*  In  addition 

to  calcium-magnesium  carbonate,  it  contains  nearly  25  percent  calcium- 
magnesium  silicate. 

Although  all  three  rock  types  were  metamorphic  in  origin,  mono¬ 
mineralic  in  appearance,  crystalline,  and  similar  in  texture  and  gram 
size,  the  shape  of  samples  of  Wausau  quartzite  were  frequently  different 
from  the  other  two  rock  types.  This  rock  comes  in  three  relatively  dis¬ 
tinct  shapes — flat,  elongated,  and  equidimensional,  whereas  anorthosite 
and  Felch  marble  were  more  generally  equidimensional  in  shape.  The 
Wausau  quartzite  shapes  are  probably  related  to  metamorphic  and  tectonic 
activities  during  the  cooling  period  of  the  Rib  Hill  mountain  region. 

Fabric  Analysis 


The  mapping  of  both  macroscopic  cracks  and  microscopic  flaws  re¬ 
vealed  that  all  three  rock  types  contained  a  large  number  of  cracks  and 
flaws.  The  distribution  was  random.  The  cracks  and  flaws  formed  a  net¬ 
work,  one  crack  or  flaw  crossing  many  others.  It  was  not  found  possible 
to  determine  the  number  or  size  of  cracks  in  any  individual  specimen. 
Figure  3  illustrates  problems  associated  with  crack  density  or  crack 
propagation  theories  of  rock  fragmentation.  This  photograph  is  typical 
of  rock  specimens  used  in  this  investigation.  Macroscopic  cracks  are 
few  in  number  and  very  little  is  known  about  their  definite  origin, 
while  microscopic  flaws  are  numerous  and  probably  associated  with  grain 
boundaries. 

It  seems  obvious  from  this  photograph  that  the  effects  of  macro¬ 
scopic  cracks  are  more  significant  than  the  effects  of  microscopic  flaws 
on  certain  rock  properties.  This  is  only  a  qualitative  conjecture 
because  it  is  not  possible  to  verify  or  quantify  this  statement  with 
experimental  evidence.  In  certain  cases,  macroscopic  cracks  are  quite 
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thin  and  their  direction  is  irregular  due  to  joining  or  intersection  of 
other  cracks.  If  there  are  microscopic  flaws  beneath  the  surface,  their 
effect  iu  unknown  The  direction  of  the  crack  propagation  is  not  related 
to  the  direction  cf  the  load.  1so,  it  is  not  possible  to  determine 
the  true  surface  area  of  an  ir  ilar  fracture  path.  Thus  it  was  con¬ 
cluded  that  it  was  not  feasible  lO  determine  the  quantitative  effect  of 
macroscopic  cracks  or  microscopic  flaws  in  this  rock  fragmentation  project 
or  to  include  such  an  effect  in  the  experimental  design. 
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BACKGROUND 


Fragment  Distributions 


A  number  of  both  empirical  and  theoretical  equations  have  been  used 
to  describe  the  size  distribution  of  broken  material.  A  brief  account  of 
the  major  equations  is  given  below.  The  data  is  usually  expressed  as  the 
cumulative  fraction  or  cumulative  percentage  by  weight  under  size,  Y, 

(o:;  over  size  1  -  Y  =  Y1)  in  relation  to  size  x. 


Three  of  the  common  equations  used  are  as  follows: 
(1)  Gates  -  Gaudin  -  Schuhmann  (8^,  6_,  _15)  : 

Y 


(-)m 

V 


where  Y  =  cumulative  fraction  finer  than  size  x, 

m  =  distribution  modulus  referring  to  the  spread  of  the  distribution, 
i.e.,  to  the  slope  of  the  line  on  log-log  paper, 

i  __  _ 2/ 


size  modulus£'  . 


—  Schuhmann  (15)  expressed  the  size  distribution  of  broken  material  by 


defining  an  extrapolated  intercept  from  the  linear  part  (fines)  of  the 
cumulative  curve.  The  extrapolation  of  the  linear  portion  of  the- curve 
to  the  100  percent  passing  ordinate  defines  the  100  percent  size  modulus;, 
k.  This  method  of  expressing  the  product  size  has  been  used  by  many 
investigators.  (Similarly,  an  80  percent  modulus  defined  as  the  same 
line  extrapolated  to  the  80  percent  passing  ordinate  has  also  been  used.) 


(2)  Rosin  -  Rammler  (14) : 


where  Y  =  cumulative  fraction  finer  than  size  x, 


b  and  n  are  similar  to  the  constants  m  and  l/km  in  equation  1, 
(3)  Gaudin  -  Meloy  (7)  ’• 


where  Y  =  cumulative  fraction  finer  than  size  x, 
xq  «*  100  percent  passing  size,  and 


r°  =  parameter  related  to  the  mean  spacing  of  the  flaws. 


; 

l 


Scale,  cm 


Gaudin  and  Meloy  suggested  that  their  equation  describes  the  distri¬ 
bution  curve  in  the  coarse  range  while  the  Gates  -  Gaudin  -  Schumann  and 
■>  Rosin  -  Rammler  equations  describe  the  distribution  curve  An  the  fine 
range.  Thus  an  overall  satisfactory  curve  fit  would  require  the  combina¬ 
tion  of  these  two  equations.  Gaudin  and  Meloy  derived  their  equation  by 
using  probability  theory. 

Bergstrom  (1A)  suggested  that  the  equation 

Y  =  [1  -  (1  -  f-)r2]qi  (4) 

o 

combines  the  fine  and  coarse  size  curve  fitting  capabilities  of  equations 
1  and  2.  This  equation  is  then  an  extension  of  equation  3.  Here 
qi  =  constant  for  improved  curve  fit, 
x  =  100  percent  passing  size,  and 

rj>  =  parameter  related  to  the  mean  spacing  of  the  flaws. 

Harris  (10)  modified  the  Bergstrom  equation  4  by  inserting  still 
another  exponent,  as  follows: 

Y  =  1  -  {1  -  [l  -  (^-)S]r3)q2  (5) 

o 

He  claimed  no  physical  basis  for  the  insertion  of  q2  in  equation  5.  He 
merely  inserted  this  term  for  the  sake  of  increased  curve  fitting  flexibility. 
It  is  supposedly  feasible  to  determine  the  values  of  r2  and  qi  or  of  s,  r3, 
and  q2  simultaneously  by  using  graphical  or  computational  methods.  However, 
these  constants  cannot  be  related  directly  to  physical  properties  of  the 
material  or  to  the  type  of  test. 

Gilvarry  (9)  derived  the  following  theoretical  equation 

Y  =  1  -  exp  [-(|)  -  (^)2  -  (f)3]  (6) 

based  on  the  premise  that  fracture  is  caused  by  stress  activated  flaws 
randomly  distributed  in  the  volume,  in  fracture  surfaces,  and  in  the  edges 
produced  by  fracture  surfaces.  The  parameters  in  equation  6  are  defined 
as  follows: 

Y  »  cumulative  fraction  finer  than  size  x, 
x  11  mean  linear  dimension  of  a  fragment, 

k  =  constant  related  to  the  mean  spacing  of  edge  flaws  with  respect 

to  x, 

j  =  constant  related  to  mean  spacing  of  surface  flaws  with  respect 
to  x,  and 

i  *»  constant  related  to  mean  spacing  of  volume  flaws  with  respect 

to  x. 


For  small  values  of  x,  the  higher  order  terms  drop  from  this  equation 
and  it  reduces  to 

Y  =  1  -  exp  [  -  ~  ]  (7) 

This  expression  is  equivalent  to  the  Rcr.in— Raamler  equation  (2)  for  n  =  1. 

When  x  becomes  small  in  equation  7 ,  the  dominant  term  of  the  series 
expansion  of  the  exponential  becomes  1  so  that 

Y  =  £  .  (8) 

k 

This  result  is  the  same  as  equation  1  for  the  case  m  =  1. 

Klimpell  and  Austin  (12)  verified  equations  1-3  by  using  statistical 
theory  to  describe  the  fragment  sizes  obtained  in  simple  fracture  of 
brittle  solids  producing  large  particles. 

Epstein  (5)  constructed  a  statistical  model  for  the  breakage  mechanism 
and  the  breakage  process  based  on  the  theory  of  probability.  He  found  that 
under  certain  conditions  it  can  be  proved  that  the  distribution  of  broken 
material  is  asymptotically  log-normal. 

From  the  literature  cited  above,  it  was  concluded  that  the  majority 
of  investigators  have  used  various  forms  or  limiting  cases  of  the  following 
functions: 


The  exponential  distribution — 

Y  =  1  -  exp  (-Xx)  , 

the  two  parameter  Weibull  distribution — 

Y  =  1  -  exp  (-x/o)a, 


the  log-normal  distribution — 


Y  -  /X  (1/B  x/2ir)  exp  [-  (log  x  -  log  ai)2/2B2]  dx,  and 


the  power  function — 


Y  =  ax 


where  Y  is  the  cumulative  fraction  finer  than  size  x,  and  where  X;  a,o; 
aj,B;  a,b  are  parameters  of  the  four  functions  respectively. 
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Two  other  distributions  used  in  statistics  were  added  to  the  list 
of  the  above  four  functions.  These  are  as  follows: 


the  two  parameter  log-Weibull  distribution — 

Y  =  1  -  exp  {-  ( sxp  x)/0)^}, 
and  the  normal  distribution — 


(13) 


Y  -  /X  (1/a  /in)  exp  {-  (x  -  Vi)2/2o2}  dx  (14) 

o 

where  3,0;  Pi  ,0  are  parameters  of  the  two  functions  respectively. 


It  was  not  feasible  to  select  any  single  equation  from  equations 
9-14  which  describes  all  the  size  distributions  of  broken  material  that 
were  found  in  this  investigation.  For  this  reason,  it  was  decided  to  use 
a  linear  regression  computer  program  incr  -porating  the  six  equations  9-14 
and  to  determine  the  coefficients  of  correlation  in  all  six  cases  (see 
DATA  ANALYSIS  section) . 


For  primary  breakage,  the  cumulative  percentage  of  weight  retained 
(or  passing)  plotted  against  sieve  size  opening  on  a  log-log  grid  or  log 
normal  grid  or  on  other  special  grids  frequently  gives  a  straight  line 
over  a  large  section  of  the  plot.  In  the  majority  of  cases  the  overall 
curve  consists  of  a  line  with  a  slope  in  the  fine  range  which  is  different 
from  and  usually  less  than  unity,  while  in  the  coarse  range  there  is  a 
curve  or  another  line  with  a  slope  usually  higher  than  unity.  In  other 
words,  the  distribution  plot  can  be  separated  into  the  plots  of  two  dis¬ 
tributions  represented  by  two  separate  equations,  one  for  the  fine  range 
and  another  for  the  coarse  range. 


There  is  no  evidence  to  support  the  hypothesis  that  the  same  fracture 
patterns  exist  throughout  all  size  ranges.  There  is  also  no  evidence  to 
support  the  hypothesis  that  the  change  in  fracture  pattern  is  exclusively 
related  to  a  change  in  the  slope  of  the  distribution  curve.  Irregularities 
encountered  in  the  slope  may  be  •'elated  to  such  things  as  plastic  deforma¬ 
tion  at  contact  points  (where  high  stress  concentration  occurs)  ,  test 
method,  specimen  geometry,  and  secondary  breakage. 


Energy  Relationships  in  Fragmentation 


As  mentioned  in  the  Fragment  Distributions  section,  the  straight  line 
representing  a  size  distribution  on  special  graph  paper  (e.g.,  a  log-log 
plot)  and  passing  through  some  of  the  data  points  does  not  go  through  all 
the  points  of  both  the  coarse  and  fine  region.  Thus  there  is  no  real 
justification  for  expressing  the  300  percent  size  modulus  k  as  an  extra¬ 
polated  point  based  on  this  line.  Consequently  there  is  no  clear  relation¬ 
ship  between  this  modulus  and  some  other  measure  of  product  distribution. 
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Instead  of  using  a  size  modulus,  a  dimensionless  mean  product  size 
is  used  in  this  report.  The  following  definition'  was  adopted: 


u  *=  ^2.  a  mean  product  size  (dimensionless) 


n 


3/ 


where:  x2  =  E  f  ^  x^  =  mean  product  size—'  (in), 

i  »  1 


-~^The  summation  is  actually  from  i  =  1  to  ®;  however  for  the  smaller  sizes 


the  terms  are  small  and  can  be  neglected. 


x.  =  average  sieve  size  (mid-point  of  sieve  class)  (in), 
_  weight  in  sieve  class 
l  -  total  weight  of  specimen  ’ 


x^  =  average  feed  size 


4/ 


•^weight  of  specimen  —  (in) . 
/  density 


— ^The  above  definition  of  xf  is  based  on  an  assumed  cubical  shape  because 
the  specimens  tended  to  be  blocky.  Other  shapes,  such  as  spheres,  could 
also  be  used  and  would  yield  different  factors  in  front  of  th*  cube  root. 


Bergstrom  (1)  conducted  a  series  of  tests  at  loading  rates  of  .10^,  10^, 
104,  and  106  lbs/min.  He  concluded  that  the  rate  of  loading  had  little 
effect  on  the  energy  -  product  relationship  for  glass  spheres.  However,  he 
reported  an  increase  in  the  energy  required  for  fracture  as  the  rate  of 
loading  was  increased.  He  also  investigated  the  magnitude  and  effect  of 
kinetic  energy  on  secondary  breakage  of  glass  spheres  jacketed  in  steel  and 
gelatin  chambers.  He  noticed  a  significant  effect  of  kinetic  energy  on 
product  size  distribution  for  two  cases.  He  also  concluded  that  the  test 
environment,  whether  a  steel  crushing  chamber  or  gelatin  chamber,  did  not 
affect  the  average  energy  requirement  for  fracture  of  spheres  of  equal 
diameter . 


In  1867,  Rittinger  postulated  that  the  energy  required  for  size  re¬ 
duction  of  a  solid  material  would  be  proportional  to  the  new  surface  area 
created  in  fragmentation.  Because  the  new  surface  is  directly  proportional 
to  the  square  of  the  new  size  and  directly  proportional  to  the  average 
number  of  new  particles  which  in  turn  is  inversely  proportional  to  the 
cube  of  the  new  size,  mathematically  Rittinger* s  hypothesis  can  be  inter¬ 
preted  as  follows: 


E  =  K(1/x2  -  1/xi) 


(15) 
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where 


E  **  energy  input  per  unit  volume , 

K  »  constant , 

«=  final  size, 
and  xi  -  initial  size. 

This  hypothesis  is  simple  to  understand  and  several  investigators 
have  used  it  to  interpret  their  experimental  data  or  to  improve  upon  this 
hypothesis. 

In  1885  Kick  postulated  that  equivalent  amounts  of  energy  should 
result  in  equivalent  geometrical  changes  in  the  size  of  pieces  of  solids. 
This  resulted  in  the  equation 


E  =  C£n  (x^^) 


where 


E  =  energy  per  unit  volume, 

C  *»  constant, 

X2  =  final  size, 
and  xi  =  initial  size. 

In  1915,  Gates  (8)  pointed  out  that  Kick’s  hypothesis  was  at  odds  with 
rigorous  experimental  evidence.  Ever  since,  no  experimental  data  has 
been  obtained  to  dispute  this  fact  or  to  show  that  Kick  s  hypothesis  is 
generally  applicable  to  size  reduction  of  rock-like  materials. 

What  is  frequently  referred  to  as  Charles'  law  (4)  in  the  crushing 
and  grinding  literature  (and  which  includes  the  above  laws  as  special  cases) 
was  first  proposed  by  Gilliland  (16)  as  follows . 

dE  =  -C  —  <16> 


where  dE  is  increment  of  energy,  x  is  particle  size,  dx  is  increment  o 
size,  and  C  and  n  are  constants.  The  integration  of  the  above  equation 
leads  to  the  following  expression: 


E  =  JX2  [-C  dx/x11]  *  -Ci  [l/x2  n  “  1  -  1/xj 
X1  <=  Cin  (xj/x2)  for  n  =  1. 


n  -  1 


]  for  n 1, 


Thus  for  n  =  1  and  2,  equation  17  takes  the  form  of  Kick's  and  Rittinger's 
laws  of  crushing  respectively. 

Bond  (2)  proposed  that  since  neither  Kick’s  nor  Rittinger's  hypothesis 
seemed  valid  for  plant  design  work,  an  energy-size  reduction  relationship 
somewhere  between  the  two  was  more  applicable.  The  fundamental  statement 
of  Bond's  mean  index  equation  is  derived  from  equation  16  with  n  1.5. 


"C  / 


K  [l/x2  *5  -  1/xi  ,5] 
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Charles  (4)  pointed  out  that  for  a  given  cumulative  distribution 

y  =  f(x)  or  dy  =  f’(x)dx  (3 

where  y  =  percent  of  weight  less  than  size  x,  the  energy  dE  required  to 
reduce  a  weight,  dy,  from  size  x^  to  x  is 

dE  =  /  [-C  dx/xn]  dy  (- 


The  energy  E  is  then  given  by 

lc  X 

E  =  /  /  (-C  dx  /x  n)  f ’ (x) dx 


x 

o  m 

where  k  is  the  size  modulus  (described  earlier  in  this  section)  and  c  is 
a  constant.  Charles  showed  that  this  reduces  to  the  form 

E  =  A/kn  '  1 

under  various  assumptions.  This  is  similar  to  one  of  the  forms  used  later 
in  this  report. 

Oka  and  Majima  (13)  analyzed  energy  requirements  in  size  reduction  of 
irregular  specimens .  Their  equation  is  as  follows : 

El  -  Ki  (X1  -W  -  x2  ~6/l!)  (22) 

where:  El  =  total  energy  required  in  size  reduction  from  feed  size  xj  to 
product  size  xj, 

Kl  =  constant, 

B  =  constant. 

For  g  =  12  or  6  this  equation  reduces  to  the  empirical  laws  of  Bond  and 
Rittinger  respectively;  i.e.,  equation  22  is  another  form  of  Charles  law. 

Intuitively,  it  is  an  appealing  concept  that  the  amount  of  energy 
utilized  in  the  fragmentation  process  is  proportional  to  the  new  surface 
area  created  (Rittinger* s  hypothesis).  However,  this  concept  does  not 
hold  in  many  cases  due  to  the  inherent  heterogeneity  and  anisotropy  of 
rocks  and  rock-forming  minerals.  A  certain  fraction  of  this  energy  is 
utilized  in  creating  a  network  of  cracks  and  flaws  in  new  particles  w  ic 
were  not  inherent  in  the  virgin  material. 

Thus,  a  summary  of  the  whole  situation  of  "laws"  and  equations  is  that 
theoretical  equations  tend  to  be  approximate  and  essentially  no  more  ac¬ 
curate  than  empirical  equations. 
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Matrix-Vector  Description  of  Single  Event  Comminution 


In  order  to  explain  the  results  of  the  fragmentation  analysis  of  this 
impact  crushing  investigation,  it  is  necessary  to  relate  the  original  sieve 
size  of  the  specimens  to  the  distribution  of  the  various  product  sieve 
sizes  obtained.  One  convenient  method  of  doing  this  is  to  use  the  breakage 
matrix  and  the  selection  matrix.  These  relate  the  input  (feed)  to  the 
output  (product)  and  are  explained  in  more  detail  below.  Further  details 
can  be  found  in  Broadbent  and  Calcott  (3)  . 


Let  a l. a  a  be  a ’decreasing  geometric  sieve  size  sequence  with 


* '  —  ' 

geometric  ratio  k,  i.e.,  a^j  = 

kai, 

with  k  <  1. 

The  breakage  matrix,  B,  is 

defined  by 

‘hi 

b12 

b!3  ‘  * 

•  bln 

b21 

b22 

b23  ‘  ‘ 

‘  b2n 

b31 

b32 

b33  ‘ 

*  b3n 

B  =  (bi;j)  = 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

ft 

bnl 

bn2 

bn3  *  ‘ 

•  brm 

The  b  ,'s  relate  feed  fragments  of  size  j  to  product  fragments  of 
size  i,  i.e. ,  b  .  is  that  proportion  of  a  fragment  between  a,  j  and  a^ 
before  breakage  tihich  falls  between  a^_^  anc^  a^  after  breakage. 


The  upper  triangular  portion  of  B  (abo’re  the  diagonal  line)  has  all 
zero  elements  (i.e.,  b  =  0  for  j  >  i)  because  a  fragment  cannot  increase 
in  size  during  comminution. 


The  selection  matrix  S  is  also  used  in  this  formulation.  This  matrix 
expresses  the  probability  of  selecting  a  set  of  size  grade  fragments. 

Thus  S  is  defined  as  follows: 


’1 

0  s 
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0  s 
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0 
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•  • 


•  • 


•  •  • 
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where  is  the  proportion  of  the  j-th  size  grade  selected  for  breakage. 
The  input  (or  feed  or  frequency)  vector  f  is  defined  as  follows: 


.  t 


where  f  is  the  proportion  by  weight  of  the  feed  fragments  between  the 
(i-1)  and  i-th  grade  size.  The  fraction  f  is  the  last  sieve  range  measured. 
Consequently  f^+j  is  the  fraction  undersize  and 

f1+f2  +  ...  +  f. + 

The  output  (product)  vector  p  is  defined  in  a  similar  fashion  as  follows: 


where  p.^  is  the  proportion  by  weight  of  the  output  fragments  between  the 
(i-1)  and  i-th  sieve  sizes. 

Putting  these  definitions  together,  it  is  seen  thet  the  effect  of  a 
selection  and  breakage  process  on  a  feed  distribution  f  is  represented  by 

p  =  [BS  +  (I  -  S)]  ? 

where  BS  ?  is  the  product  selected  and  broken  and  (I— S)  ?  is  that  part  of 
the  feed  not  selected  for  breakage.  This  simplifies  to 

p  »  b? 

in  the  present  case  because  S  is  the  identity  matrix  I  -  all  the  specimens 
prepared  are  selected  for  comminution. 


A  further  simplification  of  the  B  matrix  is  possible  in  many  actual 
fragmentation  processes  and  it  was  necessary  to  assume  it  to  be  valid 
for  the  present  investigation.  Let  the  matrix  B  be  of  such  a  nature 
that  fragments  of  every  size  are  broken  in  the  same  way  and  the  product 
depends  only  on  a  scale  factor.  In  other  words,  "scaling"  is  valid. 
Then,  because  the  sieve  sizes  have  been  chosen  to  be  in  geometric  ratio, 
it  is  seen  that 


aj  (k)  (k  <  1}  i  =  1,  2,  . ..,  n). 


In  this  investigation  k  =  1/  ^2  • 


Then  it  follows  that 


=  bi  -  j  +  1,  1  (i  -  j) 


so  that  the  matrix  B  is  completely  determined  by  its  first  column  and  can 
be  written  in  the  new  form: 


[bj  0  0  ...  O' 


b2  bj  0  .  0 


b2  b j  .  .  •  0 


B'  =1  . 


.  with  p  =  B1  ?  (24) 


Lb  b  .  b  „  . 

n  n-1  n-2 


where  elements  along  any  given  down  diagonal  are  all  the  same. 


Writing  out.  the  matrix-vector  multiplication,  the  following  equations 
are  equivalent  to  the  more  concise  notation  of  the  matrix— vector  formula¬ 
tion  : 


»i  "  bi£i 


Vl  +  blf2 


bf.+b  ,f9  +  ...  +  b.f 

n  1  n-1  L  in 


Ate 
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The  importance  of  the  simplification  of  B  to  B'  in  this  investiga¬ 
tion  is  that  all  the  non-zero  elements  of  the  matrix  can  be  obtained  by 
simply  finding  the  elements  of  the  first  column  (or  last  row). 

In  order  to  substantiate  the  assumption  that  the  more  general 
breakage  matrix  B  can  in  fact  be  simplified  as  described  above,  it 
would  be  appropriate  co  start  with  the  n  different  sizes  of  specimens 
and  to  fragment  a  group  of  each  size  to  find  out  if  the  scaling  is  in 
fact  valid.  In  practice  it  was  found  during  the  first  series  of  drop 
tests  that  it  was  not  possible,  within  the  constraint  on  the  number  of 
tests  that  could  be  carried  out,  to  obtain  consistent  results  of  the 
different  size  fragmentation  products  when  more  than  one  initial  size 
set  of  irregular  specimens  was  used.  Disc  specimens  were  also  used  and 
it  appeared  that  more  consistent  breakage  with  fewer  specimens  could  be 
obtained . 

Thus  one  of  the  objectives  of  this  investigation  was  to  obtain  the 
first  column  of  the  breakage  matrix  (equation  24) .  The  values  of  the 
b  's  exhibited  in  appendices  A1  -  AS  are  the  values  ultimately  obtained. 
It  is  then  assumed  that  the  other  columns  of  B/  can  be  obtained  from  the 
b^'s  as  shown  in  equation  24. 


DATA  ANALYSIS 


Drop  Tests 


Several  series  of  drop  tests  were  carried  out  to  study  the  effects 
of  size  and  shape  so  as  to  provide  preliminary  data  for  the  design  of 
impact  pendulum  experiments. 

The  first  series  of  drop  tests  was  designed  to  investigate  the 
effects  of  shape  and  size  on  the  value  of  X,  the  exponential  parameter 
in  the  exponential  distribution: 

Y  -  e-Xx 

During  this  phase  of  the  investigation,  the  exponential  law  was  being 
used  as  a  base  law  for  the  analysis.  From  this  experiment  and  later 
drop  tests,  it  was  found  that  the  exponential  distribution  was  not  a 
good  model. 

As  a  basis  for  this  preliminary  investigation,  random  shaped  speci¬ 
mens  were  used  to  simulate  the  feed  vector  in  which  all  the  specimens 
have  the  same  initial  sieve  size.  Sixty-four  tests  were  conducted  on 
specimens  varying  between  500  and  1 ,0b0  gin  to  select  the  itwwber  tests 

in  the  final  test  series. 

From  sieve  analysis  results  and  calculated  X  values,  it  was  concluded 
that  a  minimum  of  30  tests  would  be  required  for  each  test  situation  to 
obtain  reproducible  results  for  estimating  a  composite  X .  This  value  of 
30  was  a  compromise  between  statistical  estimation  of  the  variability  of 
the  mean,  and  a  need  to  keep  the  total  number  of  tests  down  to  a  reasonable 
total.  Sequential  testing  was  also  considered  in  order  to  keep  the  number 
of  samples  to  a  minimum;  however  it  was  felt  that  the  average  sample 
number  in  a  sequential  test  would  have  to  be  interpreted  as  30  so  that 
no  savings  would  be  achieved. 

Data  analysis  of  the  first  64  tests  also  revealed  a  high  scatter  in 
the  value  of  the  distribution  parameter  X.  Therefore,  a  second  test 
series  was  designed  to  investigate  the  effect  of  shape  and  size  on  test 
results. 

To  investigate  the  shape  factor  in  this  second  test  series,  the 
specimens  were  divided  into  three  groups  based  on  their  geometrical  con¬ 
figuration  -  equidimensional ,  elongated,  or  flat.  To  incorporate  the  ef¬ 
fect  of  volume,  the  weight  of  specimens  in  each  category  was  varied  from 
100  to  7,000  gm.  All  specimens  were  dropped  from  a  constant  height  of 
35  ft. 

Results  of  the  experiments  indicated  a  high  scatter  in  the  indi¬ 
vidual  values  of  X.  The  data  also  indicated  that  equidimensional  speci¬ 
mens  show  a  lower  probability  of  breakage  than  do  elongated  specimens 


which,  in  turn,  show  a  lower  probability  of  breakage  than  the  flat  speci¬ 
mens.  The  larger  specimens  also  showed  higher  probability  of  breakage 
than  the  smaller  specimens  for  comparable  shapes.  These  results  are 
summarized  in  table  2. 

The  values  of  X  for  each  weight  class  also  increase  as  the  proba¬ 
bility  of  breakage  increases.  This  is  to  be  expected  because  1/A  is  the 
mean  of  the  exponential  distribution;  i.e.,  the  average  size  decreases 
as  the  probability  of  breakage  increases. 

In  the  final  series  of  drop  tests,  one  hundred  and  eighty  specimens 
(ninety  of  Wausau  quartzite  and  ninety  of  anorthosite)  were  fragmented. 
Thirty  of  each  rock  type  were  dropped  from  heights  of  25,  30,  and  35  ft. 
The  fragments  from  each  specimen  that  broke  were  then  sieved  to  obtain 
the  weight  of  fragments  of  each  specimen  in  each  sieve  class. 

To  analyze  the  sieve  data,  the  cumulative  sieve  percentages  were 
obtained  and  plotted  for  each  specimen  that  broke  -  all  plots  for  one 
drop  height  being  shown  on  the  same  graph.  It  was  evident  from  these 
three  graphs  that  the  sample-to-sample  variation  was  so  great  that  taking 
an  average  (forming  a  composite)  of  all  the  samples  (number  of  samples  - 
n  <  30)  would  yield  a  better  description  of  all  the  samples  and  conse¬ 
quently  show  better  the  general  trend  of  the  fragmentation  process. 

Consequently  the  next  step  taken  was  to  form  a  composite  distribu¬ 
tion  by  adding  the  n  individual  weights  together  for  each  sieve  class 
and  then  finding  the  percentage  and  cumulative  percentage  for  each  sieve 
class.  This  was  carried  out  for  each  of  the  three  drop  heights  (see 
Figs.  Aa  and  Ab) . 

Another  method  of  forming  a  composite  distribution  for  each  drop 
height  is  to  find  the  fraction  by  weight  for  each  individual  specimen 
for  each  sieve  class  and  to  average  these  n  fractions  for  each  sieve 
class. 

The  first  composite  distribution  described  above  is  equivalent  to 
the  distribution  that  would  be  obtained  by  breaking  all  the  original 
size  rocks  together.  This,  in  effect,  simulates  a  multiple  rock  single¬ 
event  crushing  process.  The  second  distribution  described  above  is  the 
statistical  distribution  that  it  would  be  natural  to  form.  Thus  both 
distributions  are  logical  from  different  points  of  view.  As  it  turns 
out,  both  distributions  are  very  nearly  the  same,  i.e.,  the  overa 
percentage  in  each  sieve  class  is  approximately  the  same  for  either 
distribution.  In  this  investigation  both  distributions  were  calculated, 
but  only  the  former  is  shown  in  the  various  figures. 


TABLE  2.  -  Experimental  data  from  low  impact  velocity  drop  tests 


Equidimensional 
Elongate 3 
Flat 


Specimen  size  (small) 

Shape 

Number  of 
replications 
(includes  both 
broken  and  un- 
jroken  specimens) 

X  for  average 
weight  of  800  gm 
(approx) 

Probability 
of  breakage 

Equidimens ional 

47 

0.060 

0.618 

Elongated 

49 

.114 

.835 

Flat 

37 

.174 

.975 

Specimen  size  (large) 


Number  of 
replications 

X  for  average 
weight  of  1,600  gm 
(approx) 

Probability 
of  breakage 

32 

0.054 

0.876 

27 

.089 

.965 

46 

.116 

1.000 

23 
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CUMULATIVE,  percent  reto 


Figure  4a 


CUMULATIVE,  percent  retained 


The  final  step  was  to  statistically  analyze  and  plot  these  cumula¬ 
tive  results  to  find  which  distribution  fitted  best.  Three  distributions 
were  used  -  the  power  law  distribution,  the  exponential  distribution, 
and  the  normal  distribution.  It  turned  out  that  the  normal  plots  had 
the  best  eerreietiuii  coefficient,  but  this  was  not  evident  until  the 
following  compile? cion  had  been  resolved. 

The  composite  set  of  data  points  for  a  given  drop  height  when 
plotted  on  normal  paper  yielded  a  set  of  data  points  which  resembled  a 
straight  line  segment  and  a  slowly  bending  curve  (see  Figs.  4a  and  4b). 

Consequently  it  was  decided  to  separate  the  sieve  data  into  two 
groups  -  "fine"  and  "coarse."  The  cut-off  point  was  selected  visually 
as  approximately  1.05  inches  for  the  anorthosite  and  1.25  inches  for  the 
Wausau  quartzite.  These  two  groups  were  then  recompiled  (renormalized) 
separately  and  each  group  was  plotted  separately.  Thus  where  there  were 
three  "dual"  segment  plots  before,  there  are  now  six  separate  sets  oc 
points  and  lines  for  each  rock  type.  These  are  shown  in  figures  4c,  ^d, 
and  Af.  These  hut  six  lins  (or.  tvorfial  were  fitted  by  tie 

least  squares  procedure.  The  correlation  coefficient  for  each  line, 
showing  goodness  of  fit,  is  also  shown  on  each  figure. 

Note  that  each  curve  or  straight  line  in  figures  4a-4f  was  plotted 
for  a  given  drop  height.  This  is  equivalent  to  plotting  for  a  given 
specific  crushing  energy  because  specific  crushing  energy  is  proportional 
to  drop  height: 

W'h 

Ec  ■  w?T  =  eh 

where:  E  =  specific  crushing  energy  (ft  lb/ft3), 

£  =  rock  density  (constant)  (lb/ft3), 

W  =  specimen  weight  (lb) , 
h  =  drop  height  (ft). 


Impact  Pendulum  Test 


For  each  separate  impact  pendulum  test,  calculations  were  made  of 
the  rebound  height  of  the  second  (rebound)  pendulum  by  two  methods.  The 
rebound  heights  by  the  two  methods  respectively  were  calculated  as  follows 


First  method: 


h  =  R  -  /R2  - 


where:  h  =  rebound  height  (in)  , 

R  =  length  of  pendulum  wire  (from  pivot  point  to  center  of  gravity 
of  rebound  pendulum)  (in) , 

x  =  measured  horizontal  travel  of  rebound  pendulum  (in) . 
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re  4c.  Normal  curves  for  fragments  ?  1.25  in  fitted  to  cumulative  size 
distributions  of  composite  data  from  drop  tests  of  irregular 
Wausau  quartzite  specimens. 
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Figure  4d.  Normal  curves  for  fragments *7  1.05  in  fitted  to  cumulative  size 
distributions  of  composite  data  from  drop  tests  of  irregular 
anorthosite  specimens. 
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Figure  4e.  Normal  curves  for  fragments  >  1.25  in  fitted  to  cumulative  size 
distributions  of  composite  data  from  drop  tests  of  irregular 
Wausau  quartzite  specimens. 
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Second  method: 


V 


where 


*L 

At 


(26) 


V  =  rebound  velocity  (ft/sec), 

Ay  =1/12  (ft), 

=  width  of  black  shield  attached  to  pendulum, 

At  =  time  (sec)  for  shield  to  travel  1/12  ft  (measured  by  timer). 


Calculations  from  the  photocell/timer  indicated  inconsistent  velocity 
values.  This  was  probably  caused  by  fluctuations  beyond  control  in  back¬ 
ground  lighting  which  affects  the  two  trigger  points  that  start  and  stop 
the  At  time  interval  on  the  digital  counter.  Consequently  these  calcula¬ 
tions  were  not  used  in  the  analysis. 


Calibration  of  Impact  Pendulum 


It  was  essential  to  calibrate  the  impact  pendulum  for  two  reasons, 

(1)  to  determine  the  energy  lost  in  the  system  and  (2)  to  check  the 
alignment  of  the  pistons  by  verifying  for  given  input  energies  the  re¬ 
peatability  of  the  piston  velocities.  High-speed  photography  was  used 
to  determine  the  relative  motions  of  the  first  pendulum,  first  piston, 
second  piston,  and  second  pendulum  with  no  specimen  present.  Five  high¬ 
speed  photographs  were  made  of  these  moving  parts.  It  was  found  that 
after  impact  the  latter  three  parts  remained  in  contact  with  each  other 
and  travelled  as  one  unit  for  a  distance  of  about  1  in.  After  this,  the 
second  pendulum  separated  from  the  two  pistons  because  of  the  roller 
bearing  friction  on  the  pistons.  The  first  pendulum  approximately  stopped 
after  impact.  Since  the  latter  three  parts  remained  in  contact  with  each 
other  briefly  after  impact,  it  was  assumed  that  the  velocities  of  all 
three  were  the  same. 


In  order  to  calculate  the  mechanical  energy  lost,  E  ,  the  following 
energy  balance  relationship  is  used  (see  sketch  below) :  1 


or 


Energy  in  =  Energy  out  +  Energy  lost, 

Wphp  =  1/2  (^2.  +  Hi)  v22  +  W4h4  +  E 


(27) 
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where 


Wi  =  weight  of  first  pendulum, 
hi  =  release  height, 

W2  =  weight  of  first  piston, 

W3  =  weight  of  second  piston, 

V2  =  velocity  of  first  aril  second  pistons, 

W4  =  weight  of  second  pendulum, 

h4  =  height  to  which  second  pendulum  ascends, 


Et  =  energy  lost  in  friction,  heat,  vibration,  etc.  with  no 
1  specimen  between  pistons. 


Equation  (27)  can  be  rewritten  as 


E  =  Wihi  -  W4I14  -  1/2  (r^  +  V22 

i»l  e  5 


(28) 


The  numerical  values  for  the  input  energy,  Wihi ,  for  the  calibration 

output  energy  1/2  (-2-  +  —2)  V22  +  W4h4 ,  and  for  energy  lost,  E^  ,  are  given 
^  8  S  1 

in  Table  3  (see  equation  28) . 


Calculation  of  Crushing  Energy 


High-speed  photographs  of  specimen  fragmentation  runs  indicated  that 
all  the  moving  parts  (first  pendulum,  first  piston,  specimen,  second 
piston,  and  second  pendulum)  remain  in  contact  immediately  after  impact. 
When  the  specimen  is  held  between  the  two  pistons,  the  motion  of  the  first 
pendulum  is  different  than  in  the  calibration  calculations  (where  no  speci¬ 
men  was  present)  in  that  it  moves  along  with  the  other  parts.  Thus,  to 
obtain  the  crushing  energy  utilized  in  breaking  the  specimen,  it  was  only 
essential  to  determine  the  velocity  of  the  second  pendulum  immediately 
after  impact  both  with  and  without  a  specimen  held  between  the  pistons. 


The  energy  balance  when  the  specimen  is  present  (see  sketch  below)  is 
given  by  the  following  equation: 


or 


Energy  in  =  Energy  out  +  Crushing  energy  +  Energy  lost, 
Wihi 


1/2  (^g  +  +  +  ^g‘ )  ^  +  W4h5  +  Ec  +  EL2  ’ 


(29) 
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where 


'-.Jl ,  W2,  W3,  W4 ,  and  hi  are  defined  above, 

Wg  =  weight  of  specimen, 

V3  =  velocity  of  first  and  second  pistons, 
hs  =  height  to  which  second  pendulum  ascends, 

Ec  =  crushing  energy, 

E  =  energy  lost  in  friction,  heat,  vibration,  etc.  with  specimen 
2  present. 

It  is  assumed  that 


•  Et  *  Et  • 

Li  L2 

Hence  when  equations  (27)  and  (29)  are  set  equal,  the  following  expression—^ 

— ^Note  also  that  W4h4  =  1/2  V22  and  W4h5  =  1/2  ^  V32. 

- - - 2 - § _ _ _ 

for  crushing  energy,  E^,  is  found: 

Ec  =  !/2  <rf  +  V22  -  1/2  (^  +  ^  ^  +  ^)  V32  +  W4h4  -  W4h5.  (30) 


The  input  energy  Wihi  is  the  same  during  both  calibration  runs  and  specimen 
runs  from  a  given  release  height;  consequently,  this  parameter  does  not 
affect  the  value  of  crushing  energy  Ec  and  so  does  not  appear  in  equation  30. 

A  computer  program  was  used  to  calculate  the  calibration  output  energy, 
and  crushing  energy,  and  the  specific  crushing  energy  for  each  specimen 
and  averages  of  these  for  each  release  height  (input  energy)  and  rock  type 
used  in  the  testing.  (See  DATA  ANALYSIS,  Impact  Pendulum  Test  section.) 

To  facilitate  analysis,  four  digital  computer  programs  were  used  to 
analyze  the  recorded  test  data  and  sieve  data  obtained  from  the  fragmented 
specimens. 


The  first  program  used  sieve  data  to  calculate  the  following:  cumula¬ 
tive  data  for  each  resultant  fragment  distribution,  the  composite  distribu¬ 
tion  from  a  number  of  these  individual  fragment  distributions,  and  the  means 
and  standard  deviations  of  the  individual  and  composite  distributions. 

All  the  sieve  data  from  the  270  fragmented  irregular  (-3.5  +3.0  in)  speci¬ 
mens  and  the  150  disc  (2  in  diameter  x  1  in)  specimens  of  the  three  rock 
types  were  run  using  this  program. 

The  second  program  was  used  to  calculate  the  rebound  velocity  of  the 
rebound  pendulum  from  both  a  photocell /timer  method  and  a  distance  measuring 
method.  (This  was  discussed  in  the  Impact  Pendulum  Test  Apparatus  section.) 
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The  third  computer  program  was  used  to  calculate  the  calibration 
output  energy,  the  crushing  energy,  and  the  specific  crushing  energy 
(crushing  energy/ specimen  volume)  for  each  specimen  and  averages  of 
these  for  each  release  height  (energy  input)  and  rock  type  used  in  the 
testing.  (The  equations  for  calculating  these  energies  were  developed 
in  the  Experimental  Program  section.)  This  data  is  summarised  in  ap¬ 
pendices  B1  -  B24. 

The  fourth  program  was  used  to  fit  various  functions  to  the  sieve 
data.  The  six  distribution  functions  given  in  equations  9-14  were  used. 
The  fitted  normal  distributions  are  shown  in  figures  5a,  5b,  5c,  5d,  5e, 
and  5f .  The  same  method  of  combining  individual  tests  into  a  composite 
distribution  is  used  as  was  used  to  obtain  the  composite  distribution 
for  the  drop  tests. 

The  output  from  these  computer  programs  was  quite  extensive  and 
consequently  is  not  totally  reproduced  in  this  report.  Instead,  a  summary 
of  the  most  pertinent  information  is  given  in  appendices  B1  -  B24  and 

Cl  -  C6. 

Included  in  appendices  B1  -  B24  are  specimen  weight,  mean  product 
size,  standard  deviation  of  product,  specimen  output  energy  (from  pen¬ 
dulum  output  velocity) ,  crushing  energy  (that  portion  of  the  input  energy 
actually  utilized  in  crushing) ,  specific  crushing  energy  (crushing  energy 
per  volume  of  specimen),  averages  of  the  above  parameters,  and  mean  and 
standard  deviation  of  the  composite  fragment  data.  (Composite  data  are 
obtained  by  pooling  the  weights  of  fragments  in  given  sieve  classes  from 
all  the  products  listed  in  that  particular  appendix.) 

Because  the  specimens,  particularly  the  irregular  ones,  were  of  dif¬ 
ferent  sizes  and  weights  even  though  they  were  of  the  same  sieve  size,  it 
was  necessary  to  compare  the  mean  product  size  and  the  product  standard 
deviation  with  an  original  size  as  follows:  The  original  weight  of  the 
specimen  was  converted  tc  volume  by  dividing  by  the  rock  density ,  and  the 
cube  root  of  this  volume  was  considered  to  be  the  specimen  size.  Dividing 
this  value  into  the  values  of  the  mean  and  standard  deviation  yielded  a 
dimensionless  mean  and  standard  deviation.  These  dimensionless  versions 
of  the  mean  and  standard  deviation  are  also  given  in  appendices  Bl  -  B24. 
The  dimensionless  mean  is  also  u«eJ  in  eavaral  fi#nrt«. 

Included  in  appendices  Cl  -  C6  are  summary  data  relating  to  various 
functions  that  were  fitted  to  the  cumulative  sieve  data.  The  columns 
headed  A  and  B  give  the  intercept  and  slope  respectively  of  the  fitted 
function  in  the  form  where  it  has  been  found  as  a  straight  line,  (e.g., 

R  =  exp  (Xx)  becomes  y  =  logR  =  Xx,  so  that  A  =  0  and  B  =  X) . 
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Figure  5a 


Normal  curves  fitted  to  cumulative  size  distributions  of 
composite  data  from  impact  tests  of  irregular  Wausau 
quartzite  specimens. 
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Figure  5b.  Normal  curves  fitted  to  cumulative  size  distributions  of 
composite  data  from  impact  tests  of  irregular  anorthosite 
specimens . 
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Normal  curves  fitted  to  cumulative  size  distributions 
of  composite  data  from  impact  tests  of  irregular 
Felch  marble  specimens. 
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Figure  5d.  Normal  curves  fitted  to  cumulated  size  distributions  of  composite 
data  from  impact  tests  of  disc  Wausau  quartzite  specimens. 
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The  other  two  columns  show  the  coefficient  of  correlation  and  the 
coefficient  of  determination.  The  correlation  coefficient  shows  how 
well  the  line  fits  the  transformed  data  points,  and  the  coefficient  of 
determination  shows  how  well  the  fitted  function  fits  the  original  points. 
Comparison  of  the  various  correlation  coefficients  is  usually  taken  to 
indicate  the  best  fit  of  those  functions  considered.  The  normal  func¬ 
tion  generally  yields  the  best  fit  and  the  power  curve  or  the  log  Weibull 
the  second  best  fit  in  these  appendices. 

Inasmuch  as  the  normal  was  taken  to  be  the  best  fit,  it  was  desirable 
to  make  a  check  on  the  suitability  of  the  normal,  so  a  separate  compari¬ 
son  between  the  mean  and  standard  deviation  of  the  product  as  calculated 
from  the  normal  distribution  and  as  calculated  from  the  composite  sieve 
data  (also  shown  in  appendices  B1  -  B24)  was  made.  Table  4  provides  the 
data  organized  for • comparison .  Generally  the  composite  mean  is  slightly 
greater  than  the  normal  mean,  but  the  composite  standard  deviation  is 
less  than  the  normal  standard  deviation.  Thus  it  appears  that  the  normal 
is  a  reasonable  fit  to  the  data. 


In  this  investigation,  the  first  attempt  at  relating  specific  energy 
and  mean  product  size  followed  along  the  lines  of  the  work  of  Bergstrom 
(2)  or  Charles  (4) ,  as  discussed  in  the  Fragment  Distributions  section. 
The  function  used  was  the  following  hyperbolic  (or  power  law)  function 


E/V  =  a/y  (31) 

where  0 

E  =  energy  (ft-lb) 

V  =  initial  specimen  volume  (ft3), 
a ,6  =  constants  related  to  the  breakage  process  and  the  rock  type, 
y  =  Xp/X^  =  dimensionless  mean  product  size, 

X  =  mean  product  sieve  size  (in) , 

=  initial  specimen  size  (in)  =  V^'3  (Vq  In  in3). 

The  mean  product  size  was  used  rather  than  the  100  percent  modulus 
(the  curve  fitted  value  of  the  ratio  of  the  largest  product  size  to  the 
initial  specimen  size)  because  the  100  percent  value  would  have  been 
between  0.9  and  1.0  in  all  cases,  making  it  difficult  to  distinguish 
between  differing  product  size  results  resulting  from  different  input 
energies . 

The  values  for  a  and  b  were  determined  by  a  digital  computer  least 
squares  fit  of  the  power  law.  These  values  are  summarized  in  a  brief 
table  included  in  figure  7.  Both  figures  6  and  7  show  the  power  law 
fitted  to  the  data  points,  figure  6  on  rectangular  graph  paper  and  figure 
7  on  log-log  paper.  The  lines  in  figure  7  show  excellent  fits  to  the 
data  points,  particularly  in  the  disc  specimen  cases.  Note  that  in 
figure  7, a  is  the  intercept,  i.e.,  the  value  of  E/Vq  when  y  =  1  (or 
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TABLE  4.  -  Comparison  of  mean  and  standard  deviation  as  calculated  from 
the  normal  distribution  and  from  composite  sieve  data 


Material 


Felch  Marble 
Felch  Marble 
Felch  Marble 

tlausau 

quartzite 

Wausau 

quartzite 

Wausau 

quartzite 

Anorthosite 

Anorthosite 

Anorthosite 

Felch  Marble 
Felch  Marble 
Felch  Marble 
Felch  Marble 
Felch  Marble 

Wausau 

quartzite 

Wausau 

quartzite 

Wausau 

quartzite 

Wausau 

quartzite 

Wausau 

quartzite 

Anorthosite 

Anorthosite 

Anorthosite 

Anorthosite 

Anorthosite 


Release 

height 


Shape 


irregular 

irregular 

irregular 


irregular 

irregular 

irregular 

irregular 

irregular 

irregular 


Mean  =  y 
(in.) 

Normal  I  Composite 


Standard 
Deviation  =  cr 
(in.) 

Normal  I  Composite 


2.304 

1.931 

1.522 


2.030 

1.631 

1.387 

2.647 

2.191 

1.530 


2.264 

1.958 

1.612 


2.052 

1.713 

1.476 

2.474 

2.198 

1.640 


0.939 

0.951 

0.806 


0.967 

0.910 

0.937 

1.094 

1.029 

0.872 


0.724 

0.851 

0.745 


0.835 

0.838 

0.852 

0.708 

0.811 

0.787 


disc 

disc 

disc 

disc 

disc 

0.860 

0.630 

0.472 

0.372 

0,322 

0.764 

0.666 

0.570 

0.492 

0.445 

0.573 

0.462 

0.390 

0.332 

0.307 

0.268 

0.281 

0.266 

0.239 

0.230 

disc 

0.700 

0.803 

0.488 

0.283 

disc 

0.535 

0.680 

0.447 

0.307 

disc 

0.397 

0.545 

0.414 

0.267 

disc 

0.318 

0.460 

0.319 

0.227 

disc 

0.166 

0.362 

0.327 

0.211 

disc 

0.856 

0.864 

0.543 

0.267 

disc 

0.662 

0.760 

0.472 

0.302 

disc 

0.556 

0.634 

0.463 

0.278 

disc 

0.420 

0.548 

0.421 

0.280 

disc 

0.339 

0.463 

0.309 

0.230 

CRUSHING  ENERGY,  E  per  VQ(ft  lb  per  ft3) 
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Figure  6.  Power  curve  plots  of  specific  crushing  energy  verus  mean 
fragment  size. 
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Figure  7.  Logarithmic  plots  and  fitted  parameters  for  pow 


X  “  X  ).  Also  note  that  each  point  in  figures  6  and  7  is  an  average 
fBr  a  given  input  energy  taken  from  10  points  for  the  disc  specimens 
and  30  points  for  the  irregular  specimens. 

Another  attempt  at  relating  specific  energy  and  mean  product  size 
was  somewhat  similar,  but  less  successful.  This  approach  consisted  of 
fitting  Charles'  law  (see  Background)  to  the  data  points.  Charles  law 

is  given  by 


E/V  =  K  Cl/(X2)  n_1  -  1/ (Xi)  n  *3 


(32) 


where : 


h 

Xl 


K,n 


E  =  specific  energy  (ft  ), 
volume  (ft3), 

mean  product  sieve  size  (in) , 

initial  product  sieve  size  (in)  (3.25  in  for  irregular  speci¬ 
mens  and  2.00  in  for  disc  specimens), 

calculated  constants  related  to  the  breaking  process  and  the 
rock  type. 


One  data  analysis  problem  here  is  that  there  is  no  linearizing  trans¬ 
formation  to  which  a  least  squares  fit  can  be  made.  Consequently,  trial 
and  error  calculations  were  used  to  find  a  value  of  n  that  gave  approxi¬ 
mately  constant  K  values  for  different  (E/V^,  X2)  data  points.  The  re¬ 
sults  of  these  calculations  are  summarized  xn  table  5. 

Note  that  the  power  law  function  and  the  Charles'  law  function  are 
not  too  dissimilar,  so  that  somewhat  similar  exponents  can  be  expected. 

The  main  difference  between  these  two  functions  is  the  subtractive  term 
in  the  latter  function.  As  pointed  out  in  the  Background  section,  neither 
of  these  functions  is  really  a  law,  but  merely  a  general  empirical  form 
which  may  or  may  not  fit  a  particular  breakage  process  and  rock  type. 

Discussion  of  Results 

An  interpretation  of  the  data  analysis  results  is  somewhat  complicated. 
The  rock  properties  shown  in  table  1  are  not  all  clearly  consistent  with 
each  other.  It  had  also  been  postulated  that  in  the  case  of  the  impact 
tests  perhaps  one  or  more  of  the  parameters  of  the  specific  energy  versus 
product  size  relationship  could  be  correlated  to  one  or  more  of  these  rock 
properties.  Also  it  had  been  postulated  that  the  parameters  themselves 
might  differentiate  between  rock  types  as  to  breakability.  Neither  postu¬ 
late  could  be  totally  proved  or  disproved. 

At. the  inception  of  this  project  it  was  planned  to  use  several  mono- 
mineralic  rocks  of  different  strengths.  Initially  Wausau  quartzite  and 
anorthosite  were  used  and  Felch  marble  was  later  added  as  this  rock  was 
expected  to  be  weaker  and  softer.  The  eventual  results  of  a  standard 
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suite  of  property  tests  (table  1)  indicated  a  consistent  trend  on  com¬ 
paring  compressive  strength  and  hardness  -  Wausau  quartzite  was  strongest 
and  hardest,  then  anorthosite,  and  Felch  marble  was  weakest  and  softest. 
However,  tensile  strength  did  not  follow  the  same  trend  -  the  low  value 
was  for  Wausau  quartzite,  which  seems  unlikely.  Also  the  calculated^ 
dynamic  Young's  modulus  turned  out  to  be  lower  than  the  static  Young  s 
modulus  for  Wausau  quartzite,  which  seems  unlikely.  In  retrospect, 
perhaps  the  block  from  which  the  Wausau  quartzite  property  test  specimens 
were  prepared  contained  more  flaws  than  expected.  Alternately  the  anor¬ 
thosite  or  Felch  marble  properties  might  be  inaccurate.  One  way  to  shed 
light  on  this  problem  would  be  to  carry  out  many  replications  of  these 
tests  with  other  samples  of  the  same  material.  This  was  not  possible 
Within  the  limited  time,  goals,  and  costs  of  this  project. 


The  drop  test  experiments  were  performed  on  irregular  (-3.5  +  3.0  in) 
specimens  with  the  drop  energy  level  being  the  only  factor  varied.  The 
analysis  of  the  distribution  results  of  these  experiments  showed  an  ap 
parent  significant  effect  of  the  level  of  energy  applied  on  the  size  dis¬ 
tribution  parameters. 


While  it  was  possible  to  make  these  observations  based  on  fitting 
the  drop  test  data  to  a  negative  exponential  distribution,  it  was  ap¬ 
parent  that  this  distribution  was  not  sufficiently  accurate  for  our  pur¬ 
poses.  One  reason  for  this  is  that  the  experimental  distribution  data 
were  always  bimodal  rather  than  unimodal,  as  originally  assumed.  Sub- 
sequent  analysis  of  the  data  showed  that  the  fine  material  can  be  fitted 
by  a  separate  normal  distribution  from  the  normal  distribution  describing 
the  major  portion  of  the  product. 


Certain  inconsistencies  were  observed  in  the  experimental  data. 

Neither  the  mean  size  nor  the  standard  deviation  of  the  product  vane 
in  a  regular  manner  with  the  applied  energy.  Also  the  slope  of  the  fitted 
lines  on  the  probability  plot  increased  with  drop  height  for  the  Wausau 
quartzite,  but  did  not  increase  for  the  anorthosite  (see  Figs.  4a  and  4b). 
One  cause  of  these  inconsistencies,  of  course,  is  that  not  all  the  ene^Y 
goes  into  breakage.  Another  cause  is  the  small  ;ample  sizes  that  resu  e 
from  lack  of  breakage,  particularly  in  the  case  of  the  25  foot  drop  heig  . 
Additional  specimens  were  not  run  to  replace  those  specimens  that  tai  e 
to  break,  because  the  impact  test  was  going  to  replace  the  drop  test. 


The  first  impact  pendulum  test  results  were  from  the  fragmentation 
sieve  data  analysis  (Figs.  5a,  5b,  5c,  5d,  5e,  and  5f).  On  comparing 
the  appropriate  lines  on  these  plots,  it  can  be  seen  that  there  was  some 
difference  in  the  breakability  of  the  three  rock  types,  but  that  this 
difference  was  not  extremely  great. 
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The  next  step  was  to  determine  whether  any  of  the  parameters  of  the 
specific  energy  versus  product  sire  relationship  might  be  correlated 
with  any  of  the  rock  properties  using  either  the  power  law  relationship 
or  the  Charles'  law  relationship.  From  the  brief  table  of  figure  7  for 
the  power  law  and  from  table  5  for  Charles'  law,  it  can  be  seen  that  the 
only  direct  relationship  is  for  the  a  (or  K)  values  respectively,  and 
then  only  for  the  disc  specimens,  not  the  irregular  specimens.  That  is, 
the  values  of  a  (or  K)  decrease  directly  as  compressive  strength  and 
hardness  decrease,  in  the  order  Wausau  quartzite,  anorthosite,  and  Felch 
marble.  If  b  (or  n)  were  constant  this  ordering  would  imply  decreasing 
breakability .  Elowever ,  b  (or  u)  is  increasing  at  tire  sosi*.  time  as  a 
(or  K)  is  decreasing.  Thus  no  clear  pattern  evolves. 

It  appears  that  the  b  (or  n)  values  are  more  important  in  assessing 
rock  breakability  than  the  a  (or  K)  values  for  the  higher  specific 
energies.  (Than  is,  the  higher  the  line  on  the  graph,  the  stronger  the 
rock  type.)  As  is  evident  from  figure  7,  in  both  the  irregular  and  disc 
Specimens  groups  of  lines,  some  of  the  lilies  cross  because  the  exponent 
b  (or  n)  eventually  "overrides"  the  intercept  a  (or  K) . 
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CONCLUSIONS  AND  RECOMMENDATIONS 
The  following  conclusions  can  be  drawn: 

1.  The  impact  pendulum  rock  fragmentation  device  proved  to  be  a  valuable 
tool  in  obtaining  fragment  data  and  energy-fragment  size  relationships. 

2.  The  rock  properties  tested  are  not  consistent  from  rock  type  to  rock 
type. 

3.  No  relationship  could  be  established  between  specific  energy  or 
fragment  distribution  and  number  and  size  of  cracks. 

4.  Drop  test  sieve  distributions  can  be  represented  by  composite  normal 
distributions. 

5.  Impact  pendulum  test  sieve  distributions  can  be  represented  by  normal 
distributions. 

6.  For  the  impact  tests,  specific  crushing  energy  versus  mean  product 
size  (for  a  given  initial  specimen  or  feed  size)  can  be  well  represented 
by  a  power  law  or  Charles'  law.  The  exponents  of  the  power  law  for  the 
different  rock  types  and  shapes  were  generally  larger  than  expected,  com¬ 
pared  to  earlier  experimental  results  reported  in  the  literature. 

7.  The  exponent  of  the  power  law  appears  to  be  the  controlling  factor 
for  the  higher  specific  energies,  i.e.,  for  higher  specific  energies, 
the  larger  the  exponent  the  stronger  the  rock. 

The  following  recommendations  are  made: 

1.  Some  additional  impact  pendulum  tests  at  slightly  higher  input  energy 
(and  hence  higher  crushing  energy)  levels  should  be  run  using  the  same 
rock  types.  Perhaps  up  to  twice  the  maximum  energy  used  in  the  current 
series  could  be  used.  (Excessive  energy  levels  would  eventually  change 
the  breakage  process  into  some  sort  of  grinding  process,  which  would  no 
longer  be  an  approximation  of  a  single  event  breakage.) 

2.  More  extensive  property  determinations  tests  should  be  run,  using  the 
same  three  rock  types. 

3.  Another  rock  type  could  be  run  if  a  suitably  stronger  and  harder 
monomineralic  rock  could  be  found. 

4.  Further  research  and  thought  need  to  be  devoted  to  other  possible 
interpretations  of  the  data  analysis. 
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lements,  b.,  (not  smoothed)  for  irregular 
elch  marble  specimens 


3.00 

2.50 
2.12 
1.75 

1.50 
1.25 
1.05 

0.875 
0.742 
-0.742  +  0.625 
-0.625  +  0.525 
-0.5250  +  0.4375 
-0.4375  +  0.3750 
-0.3750  +  0.3125 
-0.3125  +  0.2630 
-0.2630  +  0.2210 
-0.2210  +  0.1850 
-0.1850  +  0  1560 
-0.1560  +  0.1310 
-0.1310  +  0.1100 
-0.1100  +  0.0930 


Input  energy 
=  164  ft-lb 

Input  energy 
=  293  ft-lb 

0.1770 

0.1439 

.1826 

.1282 

.2589 

.1411 

.1730 

.1824 

.0566 

.0801 

.0449 

.0918 

.0295 

.0408 

.0166 

.0536 

j  *0081 

.0195 

.0102 

.0224 

.0102 

.0167 

.0050 

.0139 

.0042 

.0108 

.0034 

.0070 

.0031 

.0063 

.0022 

.0058 

.0016 

.0035 

.0014 

.0035 

.0012 

.0033 

.0010 

.0028 

.0009 

.0021 

.9916 

.9795 

.0372 

.0213 

.0203 

.0156 

.0113 

.0102 

.0088 

.0057 

.0052 

.0049 

.0041 

.0033 

.9684 


Table  A4  -Breakage  matrix  elements,  b^,  (not  smoothed)  for  disc-shaped 
(-2.12  +  1.75  in)  Wausau  quartzite  specimens 


Sieve 

size 

_ in _ 

-2.12  +  1.75 
-1.75  +  1.50 
-1.50  +  1.25 
-1.25  +  1.05 
-1.050  +  0.875 
r0.875  +  0.742 
-0.742  +  0.625 
-0.625  +0.525 
-0.5250  +  0.4375 
-0.4375  +  0.3750 
-0.3750  +  0.3125 
-0.3125  +  0.2630 
-0.2630  +  0.2210 
-0.2210  +  0.1850 
-0.1850  +  0.1560 
-0.1560  +  0.1310 
-0.1310  +  0.1100 
-0.1100  +  0.0930 


Input  energy 
=  164  ft-lb 

0.0000 
.0000 
.  .0000 
.0506 
.5110 
.0643 
.0489 
.0195 
.0391 
.0245 
.0217 
.0181 
.0204 
.0150 
.0194 
.0135 
.0108 
.0094 


b1  -  i  =  1,  2 . 18 

Input  energy  Input  energy  Input  energy  Input  energy 
=  223  ft-lb  =  293  ft-lb  =  365  ft-lb  =  455  ft-lb 


0.0000 

.0000 

.0000 

.0477 

.2609 

.1087 

.0752 

.0581 

.0465 

.0463 

.0422 

.0317 

.0278 

.0182 

.0194 

.0201 

.0171 

.0140 


0.0000 

.0000 

.0000 

.0000 

.1312 

.0401 

.1225 

.1072 

.0835 

.0542 

.0474 

.0296 

.0383 

.0305 

.0241 

.0281 

.0212 

.0180 


0.0000 

.0000 

.0000 

.0000 

.0373 

.0396 

.1114 

.0923 

.0896 

.0785 

.0697 

.0528 

.0505 

.0302 

.0297 

.0273 

.0253 

.0215 


0.0000 

.0000 

.0000 

.0000 

.0265 

.0104 

.0250 

.0614 

.0581 

.0612 

.0656 

.0673 

.0567 

.0454 

.0422 


Table  A5  -  Breakage  matrix  elements,  b.,  (not  smoothed)  for  disc-shaped 
(-2.12  +1.75  in)  anorthosite  specimens 


Sieve 

size 

_ in _ 

-2.12  +  1.75 
-1.75  +  1.50 
-1.50  +  1.25 
-1.25  +  1.05 
-1.050  +  0.875 
-0.875  +  0.742 
-0.742  +  0.625 
-0.625  +  0.525 
-0.5250  +  0.4375 
-0.4375  +  0.3750 
-0.3750  +  0.3125 
-0.3125  +  0.2630 
-0.2630  +  0.2210 
-0.2210  +  0.1850 
-0.1850  +  0.1560 
-0.1560  +  0.1310 
-0.1310  +  0.1100 
-0.1100  +  0.0930 


Input  energy 
=  164  ft-lb 

0.0000 

.0000 

.0000 

.0952 

.6126 

.0318 

.0200 

.0119 

.0228 

.0165 

.0209 

.0138 

.0162 

.0127 

.0115 

.0102 

.0100 

.0086 


Input  energy 
=  223  ft-lb 

0.0000 

.0000 

.0000 

.0658 

.4115 

.0777 

.0598 

.0504 

.0356 

.0277 

.0329 

.0208 

.0253 

.0155 

.0171 

.0164 

.0154 

.0124 


Input  energy 
=  293  ft-lb 

0.0000 

.0000 

.0000 

.0000 

.2289 

.0971 

.1589 

.0615 

.0696 

.0393 

.0389 

.0280 

.0340 

.0215 

.0209 

.0187 

.0188 

.0154 


Input  energy 
=  365  ft-lb 

0.0000 

.0000 

.0000 

.0000 

.1595 

.0487 

.0891 

.1217 

.0565 

.0666 

.0497 

.0440 

.0378 

.0267 

.0280 

.0281 

.0256 

.0209 


Input  energy 
=  455  ft-lb 


J 


0.0000 

.0000 

.0000 

.0000 

.0324 

.0516 

.1452 

.0732 

.0937 

.0648 

.0694 

.0610 

.0530 

.0333 

.0343 

.0296 

.0254 

.0223 


Table  A6  -  Breakage  matrix  elements,  b.,  (not  smoothed)  for  disc-shaped 
(-2.12  +1.75  in)  Felch  marble  specimens 


>  •  •  •  I 


2.12  +  1.75 
1.75  +  1.50 
1.50  +  1.25 
1.25  +  1.05 
1.050  +  0.875 
0.875  +  0.742 
0.742  +  0.625 
0.625  +  0.525 
0.5250  +  0.4375 
0.4375  +  0.3750 
0.3750  +  0.3125 
0.3125  +  0.2630 
0.2630  +  0.2210 
0.2210  +  0.1850 
0.1850  +  0.1560 
0.1560  +  0.1310 
0.1310  +  0.1100 
-0.1100  +  0.0930 


Input  energy 
=  164  ft-lb 

Input  energy 
=223  ft-lb 

Input  energy 
=  293  ft-lb 

Input  energy 
=  365  ft-lb 

o.ooco 

0.0000 

0.0000 

0.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000  . 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.4954 

.1378 

.0684 

.0938 

.1152 

.0360 

.0858 

.1256 

.1008 

.1258 

.0411 

.0737 

.1106 

.1205 

.0412 

.0528 

.0811 

.0967 

.0312 

.0565 

.0667 

.0733 

.0330 

.0445 

.0428 

.0514 

.0218 

.0297 

.0437 

.0510 

.0210 

.0297 

.0369 

.0498 

.0136 

.0204 

.0241 

.0305 

.0107 

.0175 

.0234 

.0275 

.0125 

.0158 

.0226 

.0274 

.0098 

.0137 

.0204 

.0235 

.0086 

.0120 

.0166 

.0201 

.9195 

.8874 

.8427 

.8019 

Input  energy 
=  455  ft-lb 


.0787 

.0882 

.0694 

.0697 

.0528 

.0391 

.0342 

.0323 

.0281 

.0224 
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Table  B-5.  -Mean  product  size  and  specific  crushing  energy  for 
material-AN,  shape-irregular,  from  height-B. 
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Table  B-12  -  Mean  product  size  and  specific  crushing  energy  for 
material-WQ,  shape-cylinder,  from  height-B. 
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APPENDIX  C— COMPARISON  OF  VARIOUS  DISTRIBUTION  FUNCTION  FITS 
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Table  C-l.  -  Comparison  of  various  distribution  function  fits  -  Wausa< 
quartzite,  irregular  shape 


Line :  Y  =  A  +  BX 

Correlation 

coefficient 

Coefficient 

of 

determination 

Distribution 

Intercept 

A 

Slope 

B 

WAUSAU  QUARTZITE-HEIGHT  A,  IRREGULAR  SHAPE,  30  TESTS 

Negative  exponential.. 

0.000 

-0.739 

-0.977 

0.899 

Log  Weibull . 

-2.423 

1.338 

0.968 

0.965 

Weibull . 

-0.569 

1.075 

0.983 

0.972 

Power  curve . 

-0.967 

0.843 

0.966 

0.988 

Log  normal . 

-0.034 

0.807 

0.949 

0.928 

Normal . 

-1.480 

1.067 

0.994 

0.993 

WAUSAU  QUARTZITE-HEIGHT  B,  IRREGULAR  SHAPE,  30  TESTS 

Negative  exponential.. 

0.000 

-0.620 

-0.930 

0.785 

Log  Weibull . 

-3.007 

1.492 

0.979 

0.977 

Weibull . 

-0.955 

1.176 

0.975 

0.949 

Power  curve . 

-1.274 

0.977 

0.991 

0.968 

Log  normal . 

-0.317 

0.809 

0.928 

0.900 

Normal . 

-1.790 

1.097 

0.997 

0.997 
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WAUSAU  QUARTZITE-HEIGHT  C,  IRREGULAR  SHAPE,  30  TESTS 


Negative  exponential.. 

0.000 

-0.406 

-0.907 

0.739 

Log  Weibull . 

-3.717 

1.602 

0.984 

0.978 

Weibull . 

-1.522 

1.249 

0.969 

0.881 

Power  curve . 

-1.736 

1.111 

0.983 

0.890 

Log  normal . 

-0.715 

0.755 

0.921 

0.838 

Normal . 

-2.098 

1.033 

0.997 

0.996 

\  j 

Table  C-2  -  Comparison  of  various  distribution  function  fits  -  anorthosite 
irregular  shape 


Distribution 


Line :  Y  =  A  +  BX 
Intercept  Slope 

A  B 


Correlation 

coefficient 


Coefficient 

of 

determination 


ANORTHOSITE-HEIGHT  A,  IRREGULAR  SHAPE,  30  TESTS 


ANORTHOSITE-HEIGHT  B,  IRREGULAR  SHAPE,  30  TESTS 


Negative  exponential. 

Log  Weibull . 

Weibull . 

Power  curve . 

Log  normal . . 

Normal . 


0.000  -0.358 

-0.856 

-3.798  1.542 

0.986 

-1.693  1.191 

0.962 

-1.878  1.071 

0.980 

-0.833  0.702 

0.905 

-2.128  0.971 

0.991 

ANORTHOSITE-HEIGHT  C,  IRREGULAR  SHAPE,  30  TESTS 


Negative  exponential 

Log  Weibull . 

Weibull . 

Power  curve . 

Log  normal . 

Normal . 


0.000  -0.233 

-0.824 

-4.545  1.630 

0.991 

-2.336  1.233 

0.947 

-2.456  1.152 

0.962 

-1.208  0.648 

0.885 

-2.418  0.913 

0.987 

0.669 

0.982 

0.818 

0.819 

0.778 

0.979 


0.996 

0.639 

0.623 

0.625 

0.961 


Negative  exponential. 

Log  Weibull . 

Weibull . 

Power  curve . 

Log  normal . 

Normal . 


0.000 

•2.892 

-0.849 

-1.206 

-0.222 

-1.753 


II  I  l  iHfci'  i .  — . -  .  - , 


I  li 


I'l 


b  ii 


i. 


I  ii 


if 


^  o 


I  *  If 


Table  C-3  -  Comparison  of  various  distribution  function  fits  -  Felch  marble, 
irregular  shape 


85 


Distribution 

Line:  Y  =  A  +  BX 

Correlation 

coefficient 

Coefficient 

of 

determination 

Intercept 

A 

Slope 

B 

FELCH  MARBLE-HEIGHT  A,  IRREGULAR  SHAPE, 

30  TESTS 

L _  .  ... 

Negative  exponential.. 

0.000 

-0.768 

-0.912 

0.652 

Log  Weibull . 

-3.138 

1.633 

0.983 

0.987 

Weibull . 

-0.900 

1.276 

0.970 

0.946 

Power  curve . 

-1.269 

1.038 

0.990 

0.962 

Log  normal . 

-0.232 

0.900 

0.915 

0.888 

Normal . 

-1.888 

1.240 

0.998 

0.997 

FELCH  MARBLE-HEIGHT  B,  IRREGULAR  SHAPE, 

30  TESTS 

Negative  exponential.. 

0.000 

-0.442 

-0.931 

0.777 

Log  Weibull . 

-3.561 

1.590 

0.980 

0.977 

Weibull . 

-1.375 

1.251 

0.973 

0.912 

Power  curve . 

-1.610 

1.100 

0.987 

0.925 

Log  normal . 

-0.618 

0.775 

0.930 

0.866 

Normal . 

-2.029 

1.050 

. 

0.998 

0.998 

FELCH  MARBLE-HEIGHT  C,  IRREGULAR  SHAPE, 

30  TESTS 

Negative  exponential.. 

0.000 

-0.324 

-0.852 

0.594 

Log  Weibull . 

-4.565 

1.816 

0.991 

0.990 

Weibull . 

-2.097 

1.385 

0.954 

0.739 

Power  curve . 

-2.261 

1.272 

0.970 

0.721 

Log  normal . 

-1.040 

0.758 

0.895 

0.714 

Normal . 

-2.453 

1.064 

0.994 

0.984 

Table  C-4  -  Comparison  of  various  distribution  function  fits 
quartzite,  disc  shape 


Wausau 


Line:  Y 

=  A 

+  BX 

Correlation 

coefficient 

Coefficient 

Distribution 

Intercept 

A 

Slope 

B 

of 

determination 

WAUSAU  QUARTZITE-HEIGHT  A,  DISC  SHAPE,  10  TESTS 


Negative  exponential.. 

0.000 

-4.175 

-0.998 

0.982 

Log  Weibull . 

-0.821 

2.750 

0.979 

0.972 

Weibull . 

1.368 

0.957 

0.992 

0.986 

Power  curve . 

0.131 

0.462 

0.993 

0.974 

Log  normal . . 

1.887 

1.034 

0.978 

0.951 

Normal . 

-0.510 

3.058 

0.992 

0.992 

WAUSAU  QUARTZITE-HEIGHT  1,  DISC  SHAPE,  10  TESTS 


Negative  exponential. . 

0.000 

-3.018 

-0.978 

0.886 

Log  Weibull . . 

-1.458 

3.170 

0.997 

0.996 

Weibull. . . 

1.016 

1.064 

0.974 

0.962 

Power  curve . 

0.060 

0.643 

0.996 

0.993 

Log  normal . 

1.404 

1.015 

0.944 

0.913 

Normal . 

-0.994 

3.125 

0.998 

0.997 

WAUSAU  QUARTZITE-HEIGHT  B,  DISC  SHAPE,  10  TESTS 


Negative  exponential.. 

0.000 

-2.089 

-0.990 

0.946 

Log  Weibull . 

-1.503 

2.684 

0.995 

0.993 

Weibull . 

0.595 

0.903 

0.974 

0.954 

Power  curve . . 

-0.098 

0.615 

0.993 

0.981 

Log  normal . 

0.902 

0.790 

0.951 

0.919 

Normal . 

-0.957 

2.411 

0.997 

0.995 

WAUSAU  QUARTZITE-HEIGHT  2,  DISC  SHAPE,  10  TESTS 


Negative  exponential.. 

0*000 

-1.728 

-0.914 

0.795 

Log  Weibull . 

-1.836 

2.539 

0.987 

0.975 

Weibull . 

0.285 

0.934 

0.942 

0.905 

Power  curve . 

-0.311 

0.663 

0.987 

0.936 

Log  normal . 

'  0.632 

0.789 

0.884 

0.857 

Normal . 

-1.194 

2.233 

0.964 

0.955 

WAUSAU  QUARTZITE-HEIGHT  C,  DISC  SHAPE,  10  TESTS 


Negative  exponential.. 

0.000 

-1.327 

-0.803 

0.427 

Log  Weibull . 

-2.244 

2.483 

0.941 

U  .863 

Weibull . . . 

-0.203 

0.885 

0.870 

0.70/ 

Power  curve . 

-0.653 

0.670 

0.946 

0.700 

Log  normal . 

0.209 

0.695 

0.780 

0.665 

Normal . 

-1.432 

2.046 

0.885 

0.800 
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Table  C-5  -  Comparison  of  various  distribution  function  fits  -  anorthosite, 
disc  shape 


Line:  Y  = 

A 

+  BX 

Correlation 

coefficient 

Coefficient 

Distribution 

Intercept 

A 

Slope 

B 

of 

determination 

ANORTHOSITE-HEIGHT  A,  DISC  SHAPE,  10  TESTS 


Negative  exponential.. 

0.000 

-2.942 

-0.967 

0.858 

Log  Weibull . . 

-1.600 

3.337 

0.995 

0.992 

Weibull . 

1.012 

1.125 

0.977 

0.968 

Power  curve . 

0.074 

0.703 

0.998 

0.996 

Log  normal . 

1.384 

1.051 

0.943 

0.922 

Normal . 

-1.096 

3.230 

0.996 

0.994 

ANORTHOSITE-HEIGHT  1,  DISC  SHAPE,  10  TESTS 


Negative  exponential.. 

0.000 

-1.955 

-0.995 

0.971 

Log  Weibull . 

-1.582 

2.719 

0.989 

0.985 

Weibull . 

0.564 

0.931 

0.986 

0.973 

Power  curve . 

-0.092 

0.657 

0.997 

0.992 

Log  normal . 

0.848 

0.790 

0.968 

0.946 

Normal . 

-0.996 

2.370 

0.997 

0.996 

ANORTHOSITE-HEIGHT  B,  DISC  SHAPE,  10  TESTS 


Negative  exponential.. 

0.000 

-1.395 

0.941 

Log  Weibull .  . 

-1.931 

2.711 

0.992 

Weibull . . 

0.191 

0.914 

0.976 

0.942 

Power  curve . 

-0.293 

0.708 

0.991 

0.964 

Log  normal . 

0.464 

0.707 

0.951 

0.912 

Normal . 

-1.199 

2.155 

0.996 

0.992 

ANORTHOSITE-HEIGHT  2,  u*.  "  SHAPE,  10  T"STS 


Negative  exponential.. 

0.000 

-1.351 

-0.359 

0.644 

Log  Weibull . 

-2.202 

2.595 

0.970 

0.929 

Weibull . 

-0.040 

0.949 

0."20 

0.816 

Power  curve . 

-0.513 

0.726 

0.974 

0  831 

Log  normal . 

0.321 

0.739 

0.847 

0.772 

Normal . 

-1.400 

2.115 

0.934 

0.889 

ANORTHOSITE- HEIGHT  C,  DISC  SHAPE,  10  TESTS 


Negative  exponential.. 

O.COO 

-0.990 

-0.772 

0.364 

Log  Weibull . 

-2.549 

2.414 

0.911 

G.  769 

Weibull .  . 

-0.562 

0.862 

0.845 

0.569 

Power  curve . 

-0.909 

0.695 

0.914 

0.547 

Log  normal . 

-0.097 

0.625 

0.753 

0.554 

Normal . 

-1.574 

1.839 

0.853 

0.719 

Table  C-6  -  Comparison  of  various  distribution  function  fits  -  Felch  marble, 
disc  shape 


Distribution 


Line;  Y 
Intercept 
A 


A  +  BX 

Slope 

B 


Correlation 

coefficient 


FELCH  MARBLE-HEIGHT  A,  DISC  SHAPE,  10  TESTS 
;  exponential..  0.000 

mil! .  -1*531 

1.082  1.130 


Negative  exponential 

Log  Weibull . 

Weibull . 

Power  curve . 

Log  normal . 

Normal . 


0.105 

1.463 

-1.049 


0.695 

1.068 

3.256 


FELCH  MARBLE-HEIGHT  1,  DISC  SHAPE,  10  TESTS 

:  exponential..  OTOOO  ‘ 

,ull. .  -1.674  3.242 

0.862  1.092 


Negative  exponential 

Log  Weibull . 

Weibull . 

Power  curve . 

Log  normal . 

Normal . 


0.026 

1.194 

-1.117 


0.718 

0.981 

3.005 


FELCH  MARBLE-HEIGHT  B,  DISC  SHAPE,  10  TESTS 

:  exponential. .  0.000 

,ull. .  -1.897  3.050 

0.492  1.030 


Negative  exponential 

Log  Weibull . 

Weibull . 

Power  curve. . .  . 

Log  normal . 

Normal . 


-0.123 

0.764 

-1.211 


0.758 

0.840 

2.561 


FELCH  MARBLE-HEIGHT  2,  DISC  SHAPE,  10  TESTS 

Negative  exponential..  0.000  "1’^n 

Log  Weibull .  -2.231  ‘  _ 

Weibull .  0.060  0.992 

Power  curve .  **n'a?s  0  7 18 

Normal .  "1-363  ~161 


0.992 

0.811 

0.718 

•’..161 


Negative  exponential 

Log  Weibull . 

Weibull . 

Power  curve . 

Log  normal . 

Normal . 


0.000 

-r.549 

-0.463 

-0.719 

-0.139 

-1.499 


2.641 

0.905 

0.797 

0.584 

1.742 


-0.971 

0.992 

0.983 

0.998 

0.953 

0.997 


-0.977 
0.193 
0.976 
0.996 
0  948 
0.997 


-0.981 

0.995 

0.979 

0.994 

0.953 

0.998 


-0.990 

0.992 

0.984 

0.993 

0.966 

0.998 


FELCH  MARBLE-HEIGHT  C,  DISC  SHAPE,  10  TESTS 
_ 77^  n.000  -0.704 


-0.994 
0.990 
0  987 
0.992 
0.973 
0.997 


Coefficient 

of 

determination 


0.885 

0.987 

0.980 

0.995 

0.939 

0.998 


0.885 

0.995 

0.962 

0.989 

0.919 

0.996 


0.926 

0.994 

0.961 

0.985 

0.927 

0.998 


0.957 

0.986 

0.963 

0.980 

0.937 

0.997 


0.973 

0.981 

0.960 

0.971 

0.938 

0.995 


